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The spectrum of meson and diquark excitations of dense quark matter is considered in the frame- 
work of the Nambu - Jona-Lasinio model with three types of massless quarks in the presense of a 
quark number chemical potential fj,. We investigate the effective action of meson- and diquark fields 
. - - both at sufficiently large values of /i > ~ 330 MeV, where the color-flavor locked (CFL) phase 

, is realized, and in the chirally broken phase of quark matter (/x < jic)- In the last case all nine 

. pseudoscalar mesons are Nambu - Goldstone (NG) bosons, whereas the mass of the scalar meson 

' nonet is twice the dynamical quark mass. In the chirally broken phase the pseudoscalar diquarks 

are not allowed to exist as stable particles, but the scalar diquarks might be stable only at a rather 
^ |. strong interaction in the diquark channel. In the case of the CFL phase, all NG bosons of the model 

D ' are realized as scalar and pseudoscalar diquarks. Moreover, it turns out that massive diquark exci- 

I tations are unstable for this phase. In particular, for the scalar and pseudoscalar octets of diquark 

resonances a mass value around 230 MeV was found numerically. In contrast, mesons are stable 
particles in the CFL phase. Their masses lie in the interval 400-^500 MeV for not too large values 
of > ^c- 
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I. INTRODUCTION 



It is well-known that at asymptotically high baryon densities the ground state of massless three-flavor QCD cor- 
responds to the so-called color - flavor locked (CFL) phase [H, Q- In this phase quarks of all three flavors as well as 
\^ ' three colors undergo pairing near the Fermi surface due to the attractive one-gluon exchange potential. The properties 
\^ • of different collective modes, including Nambu - Goldstone (NG) bosons, of the CFL phase were studied already in 
CN I the framework of weak-coupling QCD At intermediate baryon densities, related to compact star physics, weak- 
ly-^ . coupling expansion of QCD is not applicable, so the description of color superconductivity, including the CFL phase, 
is more adequate in the framework of effective theories for the low energy QCD region. In particular, since massless 
C — , excitations might play an important role in different transport phenomena, such as cooling processes of neutron stars 
' etc, different chiral type effective theories for the pseudoscalar NG bosons of the CFL phase are usually used (see, 
e.g., |,i). 

Another effective theory approach is based on the Nambu ~ Jona-Lasinio (NJL) models. Since any NJL model 
contains the microscopic quark degrees of freedom, it is especially convenient for the investigation of dynamical 



^ . processes in dense baryonic matter. In particular, in the three-flavor NJL model the CFL effect was already considered, 
e.g., in 0, (see also the review fi\), where some aspects of the phase structure of dense quark matter were discussed, 
including the influence of the s-quark bare mass, color- and electric charge neutrality conditions, external magnetic 
field, etc. In addition, in 0, the properties and structure of NG bosons of the CFL phase were considered in the 
framework of NJL models. 

One of the most noticeable differences between color superconductivity phenomena with three and two quark species 
is that the CFL effect is characterized by a hierarchy of energy scales. As it was established in different approaches 
quoted above, at the lowest scale lie NG bosons, which dominate in all physical processes with energy smaller, than 
the superconducting gap A. Evident contributors at higher energy scales are quark quasiparticles, which in the CFL 
phase have an energy greater than A. However, up to now we know much less about other excitations, whose energy 
and mass are of the order of A in magnitude. Among these particles are ordinary scalar and pseudoscalar mesons, 
massive diquarks etc, i.e. particles which might play an essencial role in dynamical processes of the CFL phase. In 
contrast, the properties of mesons and diquarks, surrounded by color superconducting quark matter, were already 
discussed in the framework of the two- flavor NJL model ^13, [HI [Hj 13 • 

In the present paper we are going to study just this type of excitations of the CFL ground state, i.e. mesons and 
massive diquarks, in the framework of the massless three-flavor NJL model. In our previous paper we have, in 
particular, obtained the equations for both scalar and pseudoscalar diquark masses in the CFL phase of the NJL 
model. There, our consideration was based on the effective action, which is a generating functional for one-particle 
irreducible Green functions. Now, using the same technique, we perform a numerical investigation of diquark masses, 
as well as of the masses of scalar and pseudoscalar mesons, vs the chemical potential in the CFL phase. Moreover, the 



2 



octet and singlet structure is established for massive mesons and diquarks (both scalar and pseudoscalar) in the CFL 
phase. In addition, the masses of diquarks and mesons in the chirally broken quark matter phase are also investigated. 



II. NJL MODEL AND ITS EFFECTIVE ACTION 



Let us consider the following NJL model with three massless quark flavors 
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G2 ^ \[fi-l'^TA\A'q][qil^TA\A'q^] + [fTA\A'q\[qTA\A'q^]y 



(1) 



A=2,b,7 A'=2,b,7 



In ID), /i > is the quark number chemical potential which is the same for all quark flavors, q = Ccf, q = q*C are 
charge-conjugated spinors, and C = 17^7" is the charge conjugation matrix (the symbol t denotes the transposition 
operation). The quark field q = qia is a flavor and color triplet as well as a four-component Dirac spinor, where 
i,a = 1,2,3. (Roman and Greek indices refer to flavor and color indices, respectively; spinor indices are omitted.) 
Furthermore, we use the notations T£i,Aa for Gell-Mann matrices in the flavor and color space, respectively (a = 

1, ...,8); To = ' If is proportional to the unit matrix in the flavor space. Clearly, the Lagrangian ([T]) as a whole 

is invariant under transformations from the color group SU(3)c. In addition, it is symmetric under the chiral group 
SU(3)lxSU(3)7? (chiral transformations act on the flavor indices of quark fields only) as well as under the baryon- 
number conservation group U(1)b and the axial group U(l)yii. ^ In all numerical calculations below, we use the 
following values of the model parameters (see, e.g., ref. [7:]): A = 602.3 MeV, GiA^ — 2.319 and G2 = 3Gi/4, where 
A is an ultraviolet cutoff parameter in the three-dimensional momentum space. 

The linearized version of the Lagrangian ([1]) contains collective bosonic fields aa{x),Tra{x), A''^^,{x), A^j^,{x) and 
looks like 
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(2) 



where here and in the following the summation over repeated indices a — 0,...,8 and A, A' — 2,5,7 is implied. 
Lagrangians ^ and ((2]) are equivalent which simply follows from the equations of motion for the bosonic fields 



aa{x) = -2Gi{qTaq), A\^,{x) = -2G2{q'' I TAXA'q), A 



TTaix) = ~2Gi{qi^^Taq), 



A%,{x)^^2G2{cf TAXA'q), 



AA'{x) = -2G2{qiYTAXA'q^), 
-2G2{qTAXA'q^). 



(3) 



In (2)-(3) aa{x), A^j^ji^,{x) and 7ra(x), A^^, (x) are scalar and pseudoscalar fields, correspondingly. 

Let us consider the flavor group SU(3)/=SU(3)i,+_r,, which is the diagonal subgroup of the chiral group. Then, 
all complex scalar diquark fields A^^,(a;) form an (3c, 3/) multiplet of the SU(3)cXSU(3)/ group, i.e. they are a 
color and flavor antitriplet. The same is true for complex pseudoscalar diquark fields A^^, (a;) which are also the 
components of an (3c, 3/)-multiplet of the SU(3)cXSU(3)/ group. Evidently, all diquarks A^^,(a:) have the same 
nonzero baryon charge. All the real <Ta{x) and 'Kaix) fields are color singlets. Moreover, the set of scalar aa{x) mesons 
is decomposed into a direct sum of the singlet and octet representations of the diagonal flavor group SU(3)/. The 
same decomposition into multiplets is true for the set of all pseudoscalar 'Kaix) mesons. Clearly, in this case the octet 
is constructed from three pions {tt^ and tt*'), four kaons {K^ , and K^) and the eta- meson (t^s), whereas the singlet 
(ryo) corresponds to the rj' meson. 

In our previous paper Q, using the intermediate bosonic Lagrangian Q and the Nambu-Gorkov formalism we 
have obtained in the one-fermion loop approximation the effective action S^s of the initial model ([T|). In terms of 
collective bosonic fields (jS)) it takes the following form: 



d^x 



AS AS* 

^AA'^AA' 



^AA'^AA' 



4Gi 



4G2 



-Tr 



scfxNG 



In Z, 



(4) 



In a more realistic case, the additional 't Hooft six-quark interaction term should be taken into account in order to break the axial 
U(l)yi symmetry However, in the present consideration we omit the 't Hooft term, for simplicity. 
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where Z is the 2 x 2-niatrix in the Nambu-Gorkov space, 

Z -- 



and the following notations are adopted 



£>+, -K 

-K* D- 



(5) 



(6) 



Besides of an evident trace over the two-dimensional Nambu-Gorkov (NG) matrix, the Tr-operation in ^ stands for 
the trace in spinor (s), flavor (f), color (c) as well as four-dimensional coordinate (x) spaces, correspondingly. Let us 
suppose that parity is conserved so that all pscudoscalar diquark and meson fields have zero ground state expectation 
values, i.e. (A^^, (x)) = and {-Kaix)) — 0. Furthermore, since at zero s-quark mass, = 0, only the competition 
between the chirally broken quark matter phase and the CFL one is relevant to the physics of dense QCD (see, e.g., 
[3]), we permit in the present consideration nonzero ground state expectation values only for aQ{x) and A^^(x) fields 
[A = 2,5,7). Namely, let {a^ix)) = cr, (A^^(a;)) = A, (A%4(x)) ee A*, where A = 2,5,7, but other boson fields 
from ([3]) have zero ground state expectation values. In the case A = 0, ct ^ quark matter is in the chirally broken 
phase, where the ground state is invariant under SU(3)cXSU(3)/xU(1)b. If A ^ 0, then the CFL phase is realized in 
the model, and the initial symmetry is spontaneously broken down to SU(3)i,-f ij+c- ^ Now, let us make the following 
shifts of bosonic fields in (gl): ao{x) aoix) + a, Af^(a;) Af^(cc) + A*, A%^{x) A%j^{x) + A, {A = 2,5,7), 
and other bosonic fields remain unshifted. (Obviously, the new shifted bosonic fields cro(a;), A^^(x) etc, now denote 
the (small) quantum fluctuations around the mean values cr, A etc of mesons and diquarks rather than the original 
fields ([3]).) In this case 



-K* , D- 



E , K 
K* , E* 



S , K 
K* , S* 



(7) 



where Kg, K*^ ^Ti^n'Sil, are the corresponding quantities ([6]), in which all bosonic fields are replaced by their own 
ground state expectation values, i.e. 0-0(2;) — > cr, iTaix) 0, A^^(a:) A, A^^, (x) — > etc, and 5*0 is the quark 
propagator matrix in the Nambu-Gorkov representation (its matrix elements Sij are given in Appendix B). Then, 
expanding the obtained expression into a Taylor-series up to second order of small bosonic fluctuations, we have 



5off (era , 71-a , A^^, , A^^*, ) = Sf^ + 5^^^ (cTa , TT^ , A^^, , A^^* ) 



where 
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^Tr,,/,wGln(5o) = -r!(a,A,A*) / d^x, 
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4Gi 
E , K 



Sn 



4G2 

E , K 
K* , E* 



+ 



(8) 



(9) 



(10) 



and r2(CT, A, A*) is the thermodynamic potential of the system. Notice that the term linear in meson and diquark 
fields vanishes in ^ due to the gap equations. 

The detailed investigations of the thermodynamic potential, performed in our previous paper [9] for the above 
accepted model parameter set, shows that at /i < /^c ~ 330 MeV the chirally broken quark matter phase with 
SU(3)cXSU(3)/xU(1)b symmetric ground state is realized in the model (in this case A = and M « 355 MeV, 
where M — \/2/3 cr is the dynamical quark mass). However, at /i > /ic the CFL phase of dense baryonic matter 
arises. In this phase M = and A varies with /i (see Fig. 2 in jOl). Below we suppose that the gap A is a real 
nonnegative number. 

In the following we will study the spectrum of meson/diquark excitations both in the CFL and chirally broken 
phases of the NJL model. Since particle masses are calculated by the use of corresponding Green functions, it is 
necessary to put a special attention to the effective action pO|) which is really a generating functional of the 



^ In spite of the fact that in the CFL phase the chiral symmetry is also broken, the notation "chirally broken phase" is used here and in 
the following for the phase without color superconductivity. 
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one-particle irreducible (IPI) two-point Green functions of mesons and diquarks both in the chirally broken and CFL 
phases, namely 

where X{x), Y{x) = aa{x), 7rb{x), A^^, (x), A^^*, (x). (To obtain a Green function (jlip in the chirally broken phase of 
quark matter (/i < /ic), one should use in the expression for the quark propagator Sq (see Appendix B), entering in 
([To)) . M Ri 355 MeV and A 0, whereas in the CFL phase (/i > /ic) all Green functions pT|) correspond to So with 
M = and values A 7^ presented in Fig. 2 of 9].) In the following, we shall say that in the theory there is a mixing 
between two different particles with corresponding fields X{x) and Y(x), if their IPI Green function TxY{x — y) is not 
identically equal to zero. Now, after performing in (jlOp the trace operation over the two-dimensional Nambu-Gorkov 
space, we obtain 



c(2) ^ 0(2) , <^(2) <^(2) 

cff *^mcsons ~ '^diquarks ' mixed' \-^^J 



where 



+ l^^scf. {^iiS^iiS + 25i2E*52iS + 822^' 822^'} , (13) 

/•AS AS* I AP AP* ■ 

^Suarks =-J d^x^^^^^-^^-^ ^ ^Tr,,,, {S,2K*S,2K* + 2S,,KS22K* + S2,KS2^K} , (14) 
'^i'Ld = ^Tr,,/, {5nE5i2A'* ^ S2i^SnK + Si2^'S22K* + 821X822^'} , (15) 

and 8ij are the matrix elements of the quark propagator matrix 8q defined in ([7]) (see also Appendix B). (Some 
necessary explanations concerning the trace-operation over coordinate space in the expressions (jl3p - (|15[) are given in 
Appendix A (see (|A4|)). It follows from these formulae that the effective action (fTS]) depends on the mesonic fields 
aa{x),Tri,{x) only, i.e. it is a generating functional of the two-point IPI Green functions of mesons. Furthermore, the 
effective action (jl4p is composed from diquark fields only, and the mixing between mesons and diquarks might occur 
because of the effective action (|T5|) . However, as a detailed analysis of the NJL model IJ) with three massless quarks 
shows, each Green function, containing mixing of mesons and diquarks, is proportional to MA. Thus, both in the 
chirally broken phase of quark matter, where A = 0, and in the CFL one, where M = 0, there is no mixing between 
mesons and diquarks, and so the term p5)) can safely be ignored in our present consideration. ^ 

Finally, note that because of the traces containing an odd number of 7^ matrices, there is no mixing between scalar 
and pseudoscalar particles in the framework of the NJL model (1), as required by parity conservation. 

III. MESON MASSES 

A. The case of the CFL phase 

Taking into account the remarks from the previous section, we have the following general expressions for the nonzero 
two-point IPI Green function of mesons which are valid in the chirally broken phase of quark matter as well as in the 
CFL one (other two-point mesonic IPI Green functions are zero in the model under consideration): 

r2o(2) r2o(2) 
r (x~v)- '^mcBons ( x - v) - '^mesons , . 

daa{y)daa{x) d'Kb(y)d'Kb{x) 

where a, 6 = 0, 1, 2, 8. In momentum space the zeros of the Fourier transformations of these functions are connected 
with meson masses. 



^ Note, if some of the current quark masses are nonzero, then due to M there arises a mixing between mesons and diquarks in the 
CFL phase. This effect is analogous to the mixing between the cr-meson and the scalar diquark in the color superconducting phase of a 
two-flavor NJL model with nonzero masses of u- and d-quarks [Til . Il2t . 
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Starting from ((TBI) , it is possible to obtain the expressions for the mesonic IPI Green functions (a, 5 = 0, 1, 2, 8) 

ZLrl Z 

+ 52l(z)7Va^l2(-2)7V* + 522(z)7VI^22(-^)7'r*] , (17) 
ro-b<Ti,(2) = - ^Tr^c/ [Sll{z)TbSll{~'Z)Tb + Si2iz)TbS2i{~z)n 

ZLrl Z 

+ S2l{z)nSi2hz)Tl + S22{z)TlS22{~zyb] . (18) 

In (fT7|) -(fT8 |) z = X — y and the matrix elements Sij{z) are presented in formulae (|Bl|) - (IB4ll . from which the Fourier 
transformations Sij (p) are directly seen. The corresponding Fourier transformations r^^^ra (p) and Try^a,, (p) look like 
(as an example, see the relation (|A7[) from Appendix A): 

r.„.Jp) = ^ + ^Tr,,/ J ^ [ Snip + qh'TaSn{qh'Ta + Si2ip + qh'T'j2i{qh'Ta 

+^2i(p + g)7Va5i2(g)7V* +^22(p + g)7'Tl^22(9)7'T-I] , (19) 

— 1 i f d^q — — — — 

^<y^<y^[p) ^ 1^^- ^^"^scf J j^;^[Sn{p + q)TbSii{q)n + Si2{p + q)TlS2i{q)Tb 

+521 (p + q)TbSi2iqyb + 522 (p + qybS22iq)Tl] . (20) 

The zeros of these functions determine the tt- and cr-meson dispersion laws, i.e. the relations between their energy 
and three-momenta. In the present paper, we are mainly interested in the investigation of the modification of meson 
and diquark masses in dense and cold quark matter. Since in this case a particle mass is defined as the value of its 
energy in the rest frame, p = (see, e.g., [ll|, [H, [13, (3), we put p — {po, 0, 0, 0) in the following. As a result, the 
calculation of two-point IPI Green functions is significantly simplified. In order to perform for the case of the CFL 
phase in (|19 p - ([20|) the trace operations over color and flavor spaces, we used the program of analytical calculations 
MAPLE. The trace over spinor space and the subsequent integration over go 7 has been performed by applying the 



ir spin 

technique elaborated in [HI, [l^. As a result, we have for the IPI meson Green functions with a = 1, 2, 8: 

Ta^.Jpo) = T^^^M^^^+A + B, (21) 



where 

A 



A (28{E+ + E^)[E+E^+E+E-] , i{E+^ + E^)[E+^E^ + E+ E'] 



B 



(2^)3 \ SE+EZ [pI - {E+ + E^r] SE+^E^ [pg _ {E+^ + E^)^] 
AjEj + E^^)[E+E^^ + E+E-] 

5E+E-^[pl-iE+ + E-^)^] 

d^q \ 8A^{E+ + E^) ^ 8A^{E+^+E^) 



(27r)3 \ SE+E^ [pI - [El + E^f] 3E+^E^ [pg _ (E+^ + E^)^] 

'2a) 



8A^{E++E-^) 



(22) 



(23) 



3E+E-^[pl-iE+ + E-^f] 

Moreover, we use in these formulae the notations (E^)^ = {E^Y + (-E^^a)^ = i^^Y + 4|Ap, E^ = E ± fi, 



E ~ y q^ + M"^, in which M is set equal to zero. It is clear from (|22l) - ((23)) that each of the Green functions ((2T|) 

depends on Pq. So the mass squared of the a-th scalar (or pseudoscalar) meson (a = 1,2, ...,8) is determined by a 
zero of the function Ta^a-aiPo) (or r7ra7r„(po)) in the Pq plane. Moreover, it is evident from (PT|) that all scalar mesons 
with a = 1, 2, 8 have the same mass in the CFL phase, thus forming an SU(3) octet of scalar mesons. In a similar 
way, all pseudoscalar mesons with a = 1,2, ...,8 form another massive SU(3) octet as well. The masses of scalar and 
pseudoscalar meson octets in the CFL phase are presented in Fig. 1 and 2, correspondingly. 

In contrast, the two-point Green functions for cro{x) and 7ro(x) mesons take another form. Indeed, 

r^o^o (po) = ^ + Q + 7^, r^o-o M = ^ + 2 - ^' (24) 



where 



^ / ^ f mEX + E^)[E+E^+E+E-] 4iE+^ + E-^)[E+^E-^ + E+E'] \ 
' (2^)n SE+E^ipl-iEi+E^)^] 3E+^E-M-iE+^+E-^r] /' 



(2^)3 \ iEXE-^[pl - {E+ + E^y] 3E+^E-^[pl - {E+^ + E'^Y^ 
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FIG. 1: The behaviour of the scalar meson masses vs fj, 
in the CFL phase. 



FIG. 2; The behaviour of the pseudoscalar meson masses 
vs p in the CFL phase. 



(in (|25p - (P^ the quantities etc are taken again at M — 0). Evidently, these mesons are singlets with respect to 
the SU(3) group, and their masses are presented also in Figs 1,2. It is clear from these figures that none of the mesons 
have a zero mass in the CFL phase, i.e. they are not the Nambu ~ Goldstone bosons (NG) of this phase. Moreover, 
one can see that in the CFL phase there is a singlet-octet mass splitting of pseudoscalar and scalar mesons, which 
however vanishes in the A = limit. Indeed, if the value A = is used in ((2T|) - (p6)) . then the IPI Green functions of 
the octet and singlet mesons are the same, i.e. the mass splitting is absent. 



B. The case of chirally broken quark matter phase 



Now, let us find the meson masses in the chirally broken phase of quark matter, i.e. at at /i < /ic, M w 355 MeV, 
and A = 0. In this case the calculation of the traces over color and flavor indices in (fT9|)-([20ll is greatly simplified, so 
in the rest frame, p — {pq, 0, 0, 0), the mesonic Green functions look like (a, b — 0, 1, 2, 8): 
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,(po) = (p^-4M2) 
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(27) 



(2^)3 E[pl-AE^y ^0^0 vFu; J J ^2^^3 Eipl-iE^]' 

(In obtaining ((27)) . the gap equation (see, e.g., 0) was used in order to eliminate the couphng constant Gi from 
the expressions ([T9l) -(|20 l) . Note also, the expressions ([27|) do not follow directly from (f2T |) -([26 |) in the A = limit.) 
Evidently, in the chirally broken phase the Green functions Ttt^tt^ (po) turn into zero at Pq = for all a = 0, 1, 8. 
It means that nine massless excitations, Nambu - Goldstone bosons, do exist in the pseudoscalar meson sector of 
the model in the chirally broken phase. This fact corresponds to a spontaneous symmetry breaking of the initial 
SU(3)LxSU(3)jixSU(3)cXU(l)BxU(l)A down to the subgroup SU(3)cXSU(3)/xU(1)b in the chirally broken phase. 
Moreover, it is clear from (|27p that in this phase there exists a nonet of scalar mesons with mass = 2M ~ 710 MeV. 

In summary, we have three main conclusions about mesons in the framework of the NJL model (1) with three 
massless quarks. Firstly, we see that nine pseudoscalar mesons are NG bosons only in the chirally broken phase of 
quark matter. In the CFL phase they cease to be NG bosons, since now they acquire finite masses (see Fig. 2). 
Secondly, the CFL breaking of the symmetry generates the octet-singlet mass splitting of mesons. Thirdly, in the 
CFL phase the mass splitting among scalar and pseudoscalar mesons occurs differently. Indeed, at /i > /ic the mass 
of the scalar octet mesons is larger than the mass of the scalar singlet meson (see Fig. 1), whereas for pseudoscalar 
mesons the situation is inverse (see Fig. 2). 



IV. DIQUARK MASSES 

A. The case of the CFL phase 

As it follows from the discussion in section |TT1 all nonzero two-point IPI Green function of diquarks both in the 
chirally broken and CFL phases might be determined through the relation 

r(2) 



Txy{x - y) = - 



S'S, 



diquarks 



dY{y)6X{x)' 



(28) 
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where the effective action iSj;q^3^j.]^g is given in ((T4)) . and X{x),Y{x) = A^^, (x), A^*^, (x) or X{x),Y{x) = 
A^^, (x) , A^*g, (x) . The last restriction again means that scalar and pseudoscalar diquarks do not mix in accord 
with parity conservation. 

As shown in our earlier paper for the case of the two-flavor NJL model, any two-point Green function of pseu- 
doscalar diquarks differs in the color superconducting phase from a corresponding Green function of scalar diquarks 
by a term which is proportional to M^, where M is the dynamical quark mass in this phase. The same is true for 
the CFL phase of the NJL model (1). So one can conclude that in the CFL phase of the model (1), where M = 0, 
each Green function of pseudoscalar diquarks is equal to the corresponding Green function of scalar diquarks, e.g., 
r^p ,A''* ~ y) — Taj^^, Ajj*^, (a^ — y), etc. Hence, to establish the spectrum of the diquark excitations of the CFL 
phase, it is enough to study the set of scalar diquarks (the mass spectrum of the pseudoscalar diquark excitations will 
be the same in the CFL phase). 

Let us consider the two-point IPI Green functions of the scalar diquarks. A more detailed analysis of the effective 
action p4)) shows that in the CFL phase, i.e. at /j, > /ic, where M — and A 0, eighteen scalar diquarks (nine 
^AA'i^) ^^'^ A^*^, (x) fields) may be divided into four sectors: s(57, 75), s(25,52), s(27, 72) and s(257). Each 
of the sectors s{AA' ,A'A), where A ^ A', is composed of A^^, (x), A^*^, (x), A^,^(x), and A^*^(a;) diquark fields, 
whereas the sector s(257) is composed of six fields, A|2(a;), A|5(a;), A^{x), A|2(a;), A|5(a;), and Aj^{x). It turns 
out that there is a mixing between diquarks entering in the same sector, whereas fields from different sectors are not 
mixed. (The analogous situation takes place for the set of pseudoscalar diquarks, which is divided into nonmixing 
sectors p{57, 75), p(25, 52), p{27, 72) and p(257).) 



1. The case ofs{AA',A'A) sectors 



Let us first study the mass spectrum of the excitations, e.g., in the sector s(57, 75). The two-point IPI Green 
functions of scalar diquarks from this sector can be obtained from ([T4)) by the relation (|28|) . In the rest frame, where 
P — (Poj Oj 0, 0), the Fourier transforms of these IPI Green functions form the following matrix [9]: 



57,75 



(Po) 



/rA|,Ai,(po) rA|,A=;(po) rA|,A?,(po) rAJ,A?-(po)\ 

rA|*A|,(po) rA-A^*(po) rA|;A?,(po) rA|-A?*(po) 

£a?,a=,(po) rA=,A-(po) rA=,A=,(po) rA=,A?-(po) 

VrA=*A|,(po) rA=*A = ;(po) rA-A?,(po) rA?-A?*(po)/ 



A C 

B C 

C A 

C B 



(29) 



where A = a + poP, B 



a 



po/3 and 



c 



(2^)3 \" 



-^2a) 



%eIe+^[pI-{eX 



+ 



4{E+r 



lOA^ 



SEX[pl-A{E+r] 



(277)3 



_fq_ 

(2^ 



(2^)3 

E+{Ei 



E 



A 

EL) 



E^, E, 



2A 



E. 



2A 

10^; 



3^^a^2a[Po ~ (-^A + E2A 

d\ 



+. ^21 



iEX[pl-A{EXY^ 



(277)3 

2A2(£;^ 



E^ 



E- 



Ej 



E 



A' 



E2A 



E. 



EIa) 



10A2 



3E+E+^[pI-{EX 



E^a)'] 



3i^+[p§-4(i?+)21 



2A r ' 



+ 



(277)3 



E^, 



E 



A' 



E. 



2A 



E. 



2A 



(30) 



(31) 



(32) 



(To obtain the above expressions for a and /3, one has to use the gap equation for A (see in 0) in order to eliminate 
the coupling constant G2 from corresponding IPI Green functions.) Evidently, in the case of mixing between particles 
the information about their masses should be extracted from the zeros of the determinant of the matrix, composed 
from corresponding IPI Green functions. So, in our case it is necessary to study the equation detr57j5(po) — 0, which 
takes the following form 



detr 



57,75 



(po) = {AB ~ C'Y = [(« - C)(a + C) - piP 



2 a2l2 



(33) 



In the p\-T^\m\e, each zero of this equation defines a mass squared of a bosonic excitation of the CFL phase ground 
state in the s(57, 75) sector. Since this sector contains four scalar diquarks, one should search for four solutions of 
the equation ([33]) in the p§-plane. Clearly, due to the structure of (|33p . this equation admits at least two different 
solutions (each being two- fold degenerate), which are given by the zeros of the expression in the square bracket. It 
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FIG. 3: The behaviour of the mass and the width of the scalar and pseudoscalar diquark octets vs fi in the CFL phase. 



was proved in [9] that a ~ C ^ Pq- Hence, the square bracket in (|33p becomes zero at the point = 0. So, in the 
s(57, 75) sector there are two massless scalar excitations, i.e. NG bosons. 

Note, the expression detT^T^^ which is in the left hand side of the equation (l33|) is a complex- valued function defined 
on the Riemann manifold composed of an infinitely large number of sheets of the variable p^. The first (physical) 
sheet is the Pq plane with the cut 4A^ < p^ along the real axis. (Just the integrals (P0)l - (|32p supply us with the values 
of the function detr57.75 on this physical sheet.) It turns out that apart from the trivial zero, pg = 0, there are no 
solutions of the equation (|33p on this sheet, so there are no stable massive diquark excitations in the s(57, 75) sector. 
Using the procedure of analytical continuation presented in [Tl| we could find for each value of the chemical potential 
fi a complex point, lying on the second sheet of pg, where the function detPsr^Ts turns into zero. Evidently, the real 
and imaginary parts of this point correspond to the mass and width of a resonance. Hence, as was pointed out above, 
in the s(57, 75) sector there appears a twicely degenerated excitation of the CFL phase, whose mass and width are 
presented in Fig. 3 as functions of p,. 

A similar situation occurs in the other four-component sectors s(25, 52) and s(27, 72). Namely, for both sectors the 
IPI Green function matrix has the form ([^5]) . Hence, in each of these sectors there are two NG bosons as well as two 
resonances with the same mass and width, given in Fig. 3. 



2. Other diquark excitations of the CFL phase 

First of all note that there are 36 (18 scalar- and 18 pseudoscalar-) diquark fields ^ in our model. So, there should 
exist at least 36 elementary diquark excitations both in the chirally broken quark matter phase and in the CFL one. 
The masses of 12 (six of them are NG bosons, the other six are massive resonances) scalar diquark excitations of the 
CFL phase were obtained in the previous section. 

As it is clear from the discussion made at the beginning of section HVAl the mass spectrum of another 12 particles, 
which are the CFL ground state excitations in the pseudoscalar diquark sectors p(57, 75), p(25, 52), and p(27, 72), is 
identical to the mass spectrum of the corresponding scalar diquarks from sectors s{AA' , A' A), where A ^ A' { see the 
previous section flV A 1[1 . Hence, in addition to the scalar excitations found in the previous section, in the CFL phase 
there exist six pseudoscalar NG bosons as well as six pseudoscalar resonances, whose mass and width are presented 
in Fig. 3. 

Concerning the CFL ground state diquark excitations from the sectors s(257) and p(257), we should note that the 
corresponding IPI Green functions form in the rest frame of the momentum representation a nontrivial 6x6 matrix 
r257(po)j which is the same both for the s(257) and p(257) diquark sectors, with a rather complicated determinant, 

detr257(po) = [{P - W){Q -Z)-{R- Tff{{2T + Rf - {2W + P){2Z + Q)}. (34) 

An exact expression for the matrix r257(po) as well as for the IPI Green functions P, Q, R, ... from (p4)) are presented 
in our previous paper where it was also shown that the equation detP257(po) = has a three- fold degenerated 
solution pg = 0. So in the diquark sectors s(257) and p(257) there are six (three scalar and three pseudoscalar) NG 
excitations, and the initial NJL model (1) as a whole has 18 NG bosons in the mass spectrum of the CFL phase. 
A more detailed consideration of the quantities P, Q, R, ... shows that P — W = A, Q — Z = B, and R — T = C, 
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where A,B,C are given in l[29)) . So, the square brackets in ([34)) are no more than detr57^75(po) presented in p3|) . ^ 
As a consequence, we see, e.g., that in the s(257) sector of the model there exist two massive resonances with the 
same mass and width, depicted in Fig. 3. Moreover, their mass and width are identical to those for the massive 
resonances from all scalar diquark sectors s{AA' , A' A). So, all these scalar diquark resonances form in total an octet 
with respect to the SU(3)L+fl-|-c group. A similar situation is valid for the pseudoscalar diquarks, where in the mass 
spectrum there is an octet of the CFL phase excitations with the same mass and width (see Fig. 3). 

Unfortunately, we did not manage to find nontrivial diquark excitations of the CFL phase from the sectors s(257) 
and p(257), corresponding to a zero of the expression F{pq) = {2T + R)^ — {2W + P){2Z + Q) that appears in the 
braces of (|34| . (Evidently, both excitations are SU(3)L+ij+c singlets.) The matter is that F{pq) is an analytical 
function on a rather complicated Riemann manifold of the variable Pq. On its first Riemann sheet there is only a 
trivial zero, = (which corresponds to the NG bosons, as it was discussed above). Due to a rather complicated 
structure of the function F{p'^), we were not able to perform its continuation onto the second Riemann sheet and get 
any information about the mass and width of the remaining SU(3)l+_r+c singlet diquark resonances from the sectors 
s(257) and p(257). 



B. The case of the chirally broken phase of quark matter 

In this phase, i.e. at < /ic, the gap A vanishes, so the matrix elements 5*12(2; — y) and 6*21(2; — y) (see (jB2p 
and (|B3p . correspondingly) of the quark propagator matrix Sq arc vanishing, too. As a consequence, the expression 

(2) 

'^diquarks '^^ ^'^'^ two-point IPI Grccn functions for diquark fields is simplified in the chirally broken phase: 



'^diquarks ~~ ^ ^ 4^ 2 "^'^ jillA >?'22A | ■ (35) 



Using this expression in ([28]) . it is possible to get for each fixed A, A' — 2, 5, 7: 

rA2,,A^^,(a;-y)= + ^^^scf {Sii{x - y)j^TAXA'S22{y - x)j^taXa'} (36) 

T^A^'A" ,{x-y)^ ~ l-TYscf {Sii{x - y)TAXA' S22{y ~ x)ta\a'} ■ (37) 

AA' AA' 4(j-2 2 

In addition, the following relations are valid: 

^^\a''^7a' - = ^i^^A'^AA. (y - ^A-^^, A- , {^-y) = Ta- ,A- (y - 2^), (38) 

and other two-point diquark IPI Green functions are identically equal to zero in the chirally broken phase. Using the 
expressions (jB5|) - (jB6P for the fermion Green functions, one can easily perform the Tr-operation over color and flavor 
indices in (P5)) - (I571) . Then, in the rest frame of the momentum representation, i.e. at p = (po; 0,0,0), we have for 
each fixed pair of A, A' = 2, 5, 7: 

rA=.,,A^,,(po) = - 16/ (^i^23]|;TW " i " ^^^^ 



Ta- ,A:;,,bo) = " ley ^-^^^-^-—-^ = — - $,(.), (40) 

where e — (po + 2p)^. Moreover, it follows from ([55]) that Fa-^^, a^*^, (po) = rAj^*^,A^^, (— Po), T^p ^^^p-^ (po) = 
F^p* ^p (—po). From the above general consideration of the diquark Green functions in the chirally broken quark 

AA' AA' 

matter phase we see that: i) scalar diquarks do not mix with pseudoscalar ones, ii) each scalar A^^, or pseudoscalar 
A^^, diquark field is mixed only with its complex conjugated one. Hence, for each pair of scalar diquarks A^^, , A^*^, 
(or pseudoscalar diquarks A^^, , A^*^, , correspondingly) we have a simple 2x2 matrix F^^/(po) of their IPI Green 
functions (it is a 2x2 matrix F^^, (po) for the system of two pseudoscalar diquarks A^^,, A^*^, , correspondingly): 



* This fact was not observed in Q , leading to an incorrect statement about the multiplet structure of the massive diquark excitations of 
the CFL phase. 
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Actually, due to the relations ([39)) - (|40p . these matrices do not depend on A, A' = 2,5,7, i.e. they are the same for 
each pair of scalar A^^, , A^*^, or pseudoscalar A^^,, A^*^, diquarks. Obviously, to obtain the diquark excitations of 

the chirally broken phase, it is sufficient to solve the equations detF^^, (po) — and detF^^, (po) — or the following 
ones 

rAlVAi,,M^4^-«i>.(.) = 0, (42) 
rA-,A^..,M^^-«i>p(e) = 0. (43) 

In the present consideration we restrict ourselves to looking only for stable diquark excitations of the chirally broken 
phase. 

Note, the functions $s(e) and $p(e) are analytical in the whole complex e-plane, except for the cut 4Af2 < e along 
the real axis. (In general, these functions are defined on complex Riemann surfaces which are to be described by 
several sheets. The integral representations for <I>s^p(e), given in ([55)) - (pn|) . define its values on the first sheet only. To 
find values of $s_p(e) on the rest of the Riemann surfaces, a special procedure of analytical continuation is needed 
(see, e.g., in [U).) Let us denote by eg and eg the solutions of the equations (|42)) and (|43)) . respectively. Of course, 
they depend on the coupling constant G2 of the diquark channel. Obviously, the stable diquark excitation corresponds 
to the root Cg which lies on the first Riemann sheet and obeys the constraint < eg < AM^. It is fulfilled only if 
H* <G2< H**, where H* and H** are defined by 



H* 



1 



4$,(4M2) 4 [AVM2 + A2 + M2 ln((A + VM^TA^)/M)] 



4$s(0) 4 [AVM2 + A2 - AP ln((A + + A^)/M)] 2 ^ ^ 

(Actually, the last equality in (|44)) . i.e. H** — 3Gi/2, is due to the gap equation for M ^ 0.) In this case eg is the mass 
squared of the stable scalar diquark in the vacuum, i.e. at fi = 0. For a rather weak interaction in the diquark channel 
(G2 < H*), eg runs onto the second Riemann sheet, and unstable scalar diquark modes (resonances) appear. Unlike 
this, a sufficiently strong interaction in the diquark channel {H** < G2) pushes eg towards the negative semi-axis of 
the first Riemann sheet, i.e. in this case eg = {MfjY < 0: where Mfj is the mass of the diquark in the vacuum. The 
latter indicates a tachyon singularity in the scalar diquark propagator, evidencing that the SU(3)l+_rxSU(3)cXU(1)_b 
symmetric ground state of the chirally broken phase is not stable (in this case there arises a deeper ground state, 
corresponding to another phase of the model, the CFL phase). A similar observation was made in the framework 
of a two-flavor NJL model, where the chirally broken quark matter phase is unstable if there is a sufficiently strong 
interaction in the diquark channel [ll|, [H, . Indeed, at a very large G2 , as it has been shown in [ITl] , the color 
symmetry is spontaneously broken even at a vanishing chemical potential. 

Let us ignore for a moment the scalar diquark sector and perform a similar analysis, based on the equation (j43p . 
for pseudoscalar diquark excitations. Then, pseudoscalar diquarks are stable excitations of the chirally broken phase 
only, if the constraint H** < G2 < H*** is fulfilled, where H** is given in (|44| and 

4$p(0) 4 [3M2A2 + A4 - 3M2AVAf 2 + A2 ln((A + ^/JVP + A^)/M)] 

(In this case the solution eg of the equation (|43)) lies inside the interval < eg < 4Af 2.) It turns out that at G2 < H** 
these excitations are unstable, whereas at H*** < G2 a tachyonic instability of the chirally broken phase appears. 

Now, combining together the above separate considerations of the scalar and pseudoscalar diquark excitations, we 
may conclude that at a rather weak interaction in the diquark channel (G2 < H*) both scalar and pseudoscalar 
diquark excitations of the chirally broken phase are resonances. If H* < G2 < H** , then, in addition to mesons, 
the scalar diquarks are stable particles in this phase (the pseudoscalar diquarks are unstable as before). Note that 
the initial massless NJL model (1) is parametrized by three independent parameters A, Gi and G2. So, one may 
expect that estimates H* and H** from (|44)) depend on A and Gi. However, as it was pointed out just after (|44p . 
the quantity H** = 1.5Gi does not depend really on the cutoff parameter A. In contrast, H* depends both on A and 
Gi. In particular, since for the parameter set accepted in sec. II we have M w 0.355 GeV, one can present in this 
case the quantity H* in the following form H* « O.66OG1. 

Having a root eg = (A/|^)2 of the equation (I42p . one can find in the case H* < G2 < H** two zeros (with respect to 
the variable pg) of the IPI Green function Fa^»^, a^^, (Po) '"^ well as four zeros of the equation detF^^, (pg) — 0. They 
provide us with the following two different mass squared of the excitations in each two scalar A^^, , A^*^, -diquark 
system: 

{M^f = (Af^ - 2^lf, {MA'f = (A/^ + 2^)2. (46) 
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In particular, for our choice of the model parameters (see sec. II) we have « 1.968M, where M « 0.355 GeV. 
Furthermore, if G2 H*^ then 2M, if G2 H*J then M'^ 0. Since there are nine scalar A^^, diquarks 

as well as nine scalar A^*^, antidiquarks in our model, we relate Ma in (|46p to the mass of the diquark nonet with 
the baryon number B = 2/3 and Ma* to the mass of the antidiquark nonet with B — —2/3. The difference between 
diquark and antidiquark masses in (|46[) is explained by the absence of a charge conjugation symmetry in the presence 
of a chemical potential ji. 

Finally, if H** < G2, then a SU(3)l+_rxSU(3)cXU(1)s symmetric ground state, i.e. the chirally broken phase of 
quark matter, is not allowed to exist in the model. The matter is that in this case tachyon singularities of the scalar 
diquark propagator alone (at H** < G2 < H***), or both of the scalar and pseudoscalar diquark propagators (at 
H*** < G2) appear. As a result, in this case the ground state of the CFL phase is always deeper in comparision with 
the ground state of the chirally broken phase. So, only the CFL phase may be realized in the model at sufficiently high 
values of the coupling constant G2 and arbitrary values of /i. As a consequence, one must expect that at G2 — *■ H** 
the critical value /ic of the chemical potential tends to zero. The fact that at G2 the diquark mass M|J, tends 

to zero may be considered as a precursor, which appears in the chirally broken phase, of the spontaneous breaking of 
the SU(3)c symmetry, taking part at G2 — H** . 

V. SUMMARY AND DISCUSSION 

In the present paper we have continued the investigation, started in our previous paper of the bosonic excitations 
(mesons and diquarks) of the dense quark matter, composed of u, d, and s quarks, at zero temperature. The consider- 
ation is performed in the framework of the massless NJL model (1), omitting the 't Hooft six-quark interaction term, 
for simphcity. In this case, the initial symmetry group of the model, i.e. SU(3)LxSU(3)flxSU(3)cXU(l)BxU(l)A 
does contain the axial U(1)a subgroup. As a result, we have shown for the model parameter set accepted in 
sec. II that at sufficiently low values oi fi, ^ < ^ 330 MeV, the chirally broken quark matter phase with 
SU(3)i+ii;xSU(3)cXU(l)B-ground state symmetry is realized and nine massless pseudoscalar mesons (which are the 
NG bosons), tt^, tt", K^, K^, K^, r]s and 77', appear. (In massless QCD, where U(1)a is broken on the quantum 
level, or in NJL models with 't Hooft interaction the '-meson is not a NG boson.) 

At ^ > /ic the original symmetry is spontaneously broken down to SU(3)L+fl,+c, and the CFL phase does occur. 
In this case, in accordance with the Goldstone theorem, eighteen NG bosons must appear in the mass spectrum of 
the model (1). (In contrast, due to the absence of the unphysical U(l)^ symmetry, only seventeen NG bosons must 
appear in massless QCD.) Considering IPI Green functions, we have found nine NG bosons in the sector of scalar 
diquark excitations. Eight of them have to be considered as non-physical, since in real QCD they supply masses to 
gluons by the Anderson - Higgs mechanism. The remaining scalar NG boson corresponds to a spontaneous breaking 
of the baryon U(1)b symmetry. The other nine NG bosons are no more pseudoscalar mesons, but now the massless 
excitations in the pseudoscalar diquark sector of the model. All that, i.e. the NG boson structure of the model (1), 
is the main result of the paper [9| which is also confirmed in the present consideration. 

Besides NG diquarks, we have proved in the existence of massive diquark excitations in the CFL phase. However, 
a detailed numerical analysis of the diquark masses vs chemical potential /i was not done there. In the present paper we 
have argued that all massive diquark excitations of the CFL phase are resonances, since the corresponding singularities 
of their Green functions in momentum space lie on the second energy Riemann sheet. Moreover, they form scalar- and 
pseudoscalar SU(3)i,-fij+c octets and singlets. The mass and width of scalar- and pseudoscalar diquark resonances 
from octets vs /x have been obtained numerically (see Fig. 3). They are avaluated around 230 MeV and 50 MeV, 
correspondingly, i.e. this quantities are at least five times smaller, than the mass and width of the scalar diquark 
resonance in the color superconducting quark matter composed of u and d quarks [ll|, [H, [13 . (Due to numerical 
difficulties, we were not able to evaluate the parameters of the above mentioned massive scalar- and pseudoscalar 
singlet diquark resonances of the CFL phase, however, we guess that their masses are of the same order in magnitude 
as the masses of the diquark octets.) 

To get a more complete view about the diquark properties in the framework of the NJL model (1), we have 
considered their masses in the chirally broken phase of quark matter, too. It follows from our analysis that i) at 
sufficiently strong interaction in the diquark channel, i.e. at G2 > H** ~ 1.56*1, the existence of this phase is 
prohibited in the framework of the NJL model (1), ii) depending on the coupling constant G2, scalar and pseudoscalar 
diquarks have different properties in this phase. Indeed, at G2 < H* , where H* is given in both types of 

diquarks are resonances. However, at H* < G2 < H** the pseudoscalar diquarks remain to be resonances, whereas 
scalar diquarks are yet stable particles. As this takes place, there is a splitting between the scalar diquark and scalar 
antidiquark masses (see (j46p ). which is explained by the violation of the charge conjugation symmetry in the presence 
of a chemical potential. (Of course, in the chirally broken phase all observable particles are colorless, so one should 
expect that colored diquarks are confined within baryons (see e.g. [13 )■ Thus, one may look at our investigation 
of the diquark masses in the chirally broken phase as an indication of the existence of rather strong quark-quark 
correlations inside baryons, which might help in a better understanding of baryon dynamics.) 

Finally, we have considered in the model (1) the masses of mesons which are stable particles in both phases. In the 
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chirally broken phase, i.e. at < Hc, all nine pseudoscalar mesons are NG bosons, whereas the nine scalar mesons 
have equal mass = 2M ~ 710 MeV (for the parameter set of the model accepted in sec. II). In the CFL phase, i.e. 
at ^ > /ic, these nonet representations of mesons, reducible with respect to the SU(3)l+_r+c group, are decomposed 
into the octet and singlet representations, each with its own mass. The reason for this octet-singlet mass splitting of 
mesons is just the color-flavor locked symmetry breaking taking place at /i > ^c- In the CFL phase the masses of 
both types of mesons vary in the interval 300-^900 MeV, when /i varies from 330 MeV to 500 MeV (see Figs 1,2). 
However, the mass splitting among the scalar and pseudoscalar mesons occurs in different ways. Indeed, as it is easily 
seen from Fig. 1, the mass of the scalar octet mesons is larger than the mass of the scalar singlet meson, whereas for 
pseudoscalar mesons the opposite situation takes place (see Fig. 2). 
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APPENDIX A: SOME FORMULAE 

The present Appendix contains some useful formulae employed in the text, 
i) Determinant: 

det ( ^ ' ^ j = det[-CB + CAC-^D] = det[DA - DBD-^C]. (Al) 



ii) Inverse matrix: 



A , B\ ^ _ ( C-^DL , -N \ _ f L ,_ -A~^BN 
C , D I " [ -L , B-^AN ) ^ \ -D-^CL , N 



(A2) 



where 

L = [AC-^ D - B]-^ , N =[DB-^A-C]-^ , L = [A- BD-^C]-^ , N = [D - CA'^ B]-^ . (A3) 

iii) Variational derivatives: Let A^ B are some operators in the coordinate space with matrix elements 
A{x^ y) =A{x — y) and B{x, y) = B{x — y), respectively. Moreover, let cr{x) and 0(a;) are some fields. Then, 



Ti[{AaB(f)} = J dxdydzduA{x, z)a{z)6{z — y)B{y,u)(/){u)S{u — x) = J dxdyA{x,y)a(i))B{y,x)4>{x). (A-4) 
It follows from fXi)) that 

r(x -y)= ^^^J^y^^^^!^^ = A{x, y)B{y, x) = A{x - y)B{y - x). (A5) 

iv) Fourier transformations: For arbitrary function F{z) it is possible to define the Fourier transformation F{p) 
by the relation 

T{p) = j d^zF{z)eP^ i. e. F{z) = j 0^Fip)e-^P\ (A6) 
Taking these relations into account, one obtains from (jA5p that 

r(p) = / ^A{q + p)Biq), (A7) 
where A{q),B(q) are Fourier transformations of the functions A{x) and B{x), correspondingly. 

APPENDIX B: QUARK PROPAGATOR MATRIX 

In the Nambu-Gorkov representation the inverse quark propagator matrix Sq^ is given in ([7|). Using the techniques, 
elaborated in [l^, [O, [H, [ill, [l3| , it is possible to obtain the following expressions for the matrix elements of the quark 
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propagator matrix 5*0 = 



Sll, Si2 

S21, S22 



(27r)4 Uo-(^ba) 9o-(^i3A) 



5i2(a: - 2/) = -lAB I ^.~^i(--v) I , / , ,, 7^ A- + — -1— I , (B2) 



where M = A± = i(l ± XilfcM). Moreover, {E^^f = + jApB^, E^^E±^i, E=^q^ + A'P and 

S = X]a=2 5 t'^'aA^. (In these and other similar expressions, qo is a shorthand notation for qq + ie ■ sgn(go), where 
the hmit e ~^ 0+ must be taken at the end of all calculations. This prescription correctly implements the role of /i as 
the chemical potential and preserves the causality of the theory.) It is clear from (|Bip - (|B4l) that all color- and flavor 
dependences in the matrix elements ^n, 5*12, S'21 and 6*22 arise only due to the matrix B. It is a 9 x 9 matrix in the 
nine-dimensional space c x / which is the direct production of color and flavor spaces. Note, in the chirally broken 
quark matter phase, where A — 0, M 0, the expressions for the matrix elements (jBl[) - (|B4[) have a simpler form. 
Namely, it is clear that in this phase 6*12 (x — y) = S2i{x — y) — and 

5n(.-.)^lcxl,x /7|^e-(-)|^ + ^j, (B5) 
J (27r)4 lqo + E+ qo - E J 



S22i. - = Ic X 1, X / (^e--^^-^ 1^ + ^ } , (B6) 



i.e. the matrix elements (jB5[) - (IB6p are proportional to the unit matrices both in the color and flavor spaces. 
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